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ABSTRACT 


A digital computer program for the geometrically nonlin- 
ear static and dynamic response of arbitrarily loaded shells 
of revolution known as SATANS is modified, removing a re- 
guirement that the applied loads be symmetric about a datum 
meridional plane. "Ihe ability to include arbitrary 
naturally-occurring initial geometric imperfections of any 
shape, location or extent is added to the program. Two 
specific axially loaded circular cylinders are analyzed 
uSing actual initial imperfections. Detailed numerical 
studies are conducted to determine the buckling loads of the 
perfect and imperfect cylinders. Several sets of boundary 
conditions are considered. An analysis procedure for imper- 
fect shells is proposed. The results indicate that a means 
to conduct a Simeseterul analysis of an arbitrarily imperfect 
shell of revolution has been developed. 

A study of the buckling behavior of the clamped shallow 
spherical cap modelled under various boundary conditions at 


the pole is included as an addendum. 
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E. ENTRODUCTION 


shells are used asvstructural elements in nearly every 
space vehicle, aircraft or other structure where light 
weight is of importance. Consequently, a large number of 
digital computer programs for the structural analysis of 
this shells have been developed. Many of these programs 
were briefly described by Bushnell in Ref. 1. An assessment 
concerning the analysis capability based on computational 
systems that have application to a broad class of shell 
structures was presented by Hartung in Ref. 2. In his as- 
sessment, Hartung discusses the limitations, scope, accuracy 
and efficiency of virtually every major shell structural 
analysis program currently available and under development. 
Of particular interest here, is the computer program SATANS 
SSeotcatic And Transient Analysis, Nonlinear, Shells =- de- 
veloped by Ball for the geometrically nonlinear analysis of 
arbitrarily loaded shells of revolution /Refs. 3 and 47. 

Birci, an section, 11, a brief description of) SATANS is 
presented including the basic theory, the program capabili- 
ties and limitations, the solution method and procedure and 
some general comments. 

Second, an wWectionwLit SATANS-—F \asnconstructed.:: Modi- 
fications to SATANS making the program more user oriented 
are discussed. Included in these modifications are a reor- 
ganization of most of SATANS into a "storable package,' 


incorporation of an equipment-independent plotting package 
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which uses a high-speed line printer, modification of itera- 
tion control, discussion of extrapolation technique and 
alteration of boundary load application procedures. 

Third, and most “important, the subject of imperfections 
is treated in section IV. Thin isotropic shells of revolu- 
tion have consistently displayed experimental buckling loads 
much lower than the theoretically predicted buckling loads. 
This disagreement has generally been attributed to (1), the 
Mathematical inability to correctly model the boundary 
conditions as they exist on the test specimen, and (2), a 
high degree of imperfection sensitivity, i.e., naturally- 
occurring imperfections in the initial shape of the shell 
cause a premature triggering of critical buckling modes and 
early shell collapse. An attempt by Dantone /Ref. 57, EO 
modify SATANS and to include initial imperfections in the 
analysis of cylindrical shells was severely hampered by a 
restriction on the use of SATANS, namely, all applied 
loading and initial conditions -- and hence, initial imper- 
fections -- must be symmetric about some meridional plane. 
Discussion of the removal of this restriction and introduc- 
tion of initial imperfections into the program are covered 
on this’ section. 

Fourth, and also in section IV, the implementation of 
the initial imperfection capability is presented. The 
mesulting program is called SATANS-II. Steps for actual 
Construction of SATANS-II from SATANS-I are presented as 


ame, Instructions’ for its’ use. 
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Fifth, in section V, actual imperfect shell analyses 
are performed on two specific axially loaded circular cylin- 
ders. The actual topography of the cylinders has been 
measured by Arbocz and Babcock /Ref. 6/7 and this data pro- 
vides the opportunity to include actual naturally-occurring 
initial imperfections into the analysis by SATANS-II. 
Detailed numerical studies are performed to determine the 
classical and consistent bifurcation buckling loads of the 
perfect cylinders and the results are compared with previ- 
ously published results. Stability solutions for the 
imperfect cylinders under six boundary conditions are 
obtained and compared with experimental values. An analysis 
procedure for imperfect shells is proposed based upon these 
results. 

Lastly, as an addendum to this report, the clamped 
shallow spherical cap is investigated under various boundary 
Conditions at the pole. While on an industrial tour with 
the Structural Analysis Division of the Palo Alto Research 
Laboratory of the Lockheed Missiles and Space Company, 
Sunnyvale, California, Ryan noted that few shell structural 
analysis programs have the capability of analyzing a shell 
of revolution which has an initial or a final pole -- 
continuous structure through the axis of revolution. This 
investigation attempted to determine the effects on the 
solution caused by substitution of a small hole or a small 


migid plug for the pole. 
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II. DESCRIPTION OF SATANS 


The computer program known as SATANS is capable of per- 
forming static and transient analyses of thin, linearly- 
elastic shells of revolution subjected to arbitrary loads 
and temperature distributions. It is based on Sanders' 
nonlinear field equations /Ref. 7/7 for the conditions of 
small strains and moderately small rotations. The program 
can be used under the following conditions: 

1. The geometric and material properties of the shell 
must be axisymmetric, but may vary along the shell meridian. 

2. The shell material must be isotropic, but the 
modulus of elasticity may vary through the thickness. 

3. Poisson's ratio must be constant. 

4, The boundaries of the shell may be closed, fixed, 
elastically restrained or free, and loaded or unloaded. 

5. The applied forces, pressures, temperatures and 
initial conditions must be symmetric about a datum meridional 


plane. 


A. GEOMETRY AND COORDINATE SYSTEM 

The coordinate system used to locate points within the 
shell is the s, 8, ft orthogonal system, where s is the 
meridional distance along a reference surface, § is the 
circumferential angle measured about the axis of revolution 
and [ is the perpendicular distance from the reference 


surface. The shell geometry and coordinate system are shown 
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in Figure 1. Positive directions for displacements, rota- 


tions, forces, moments and loads are shown in Figure 2. 


B. METHOD OF SOLUTION 

The goveuming partial differential equations are reduced 
to an infinite number of sets of four second order differen- 
tial equations in the meridional and time coordinates by 
expanding all dependent variables in a sine or cosine series 
in terms of nO, where n is the Fourier index (harmonic) and 
6 is the circumferential coordinate. The sets are uncoupled 
by utilizing appropriate trigonometric identities and by 
treating the nonlinear coupling terms as ‘pseudo loads,' 
The meridional derivatives are replaced by conventional 
central finite difference approximations, and the displace- 
ment accelerations are approximated by the implicit Houbolt 
backward differencing scheme /Ref. 87 « This leads to sets 
of algebraic equations in terms of the dependent variables 
and the Fourier index or harmonic. At each load or time 
step, an estimate of the solution is obtained by extrapola- 
tion from the solutions at previous load or time steps. The 
sets of algebraic equations are repeatedly solved using 
Potters' form of Gaussian elimination /Ref. a. The pseudo 
loads are recomputed after each solution and the procedure 
continues until the solution converges. 

In the static analysis, a load-displacement history is 
obtained by an automatic incrementing of applied loadings 
in equal steps until the number of iterations required to 


achieve convergence in any harmonic exceeds a prescribed 
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limit. The load increment is then reduced by a factor of 

five and the procedure continues until a prescribed number 
of load steps, or a prescribed number of load step reduc- 

tions have been made. 

Since the method of solution is based on the pseudo load 
approach, the shell reacts equally, ina linear fashion, to 
any change in either the applied load or the pseudo load. 
Thus, failure of the solution to converge in any harmonic 
can be attributed to two types of nonlinear behavior, both 
Growhich are illustrated in Figure 3. Referring to. the 
figure, curve A experiences an inflection point or a maximum. 
Either condition causes a zero or nearly zero slope and a 
convergence failure. This curve is typical of softening 
fenlinearity. Curve B, on the other hand, is typical of a 
stiffening nonlinearity and a convergence failure will occur 
when the slope becomes large. It is necessary, then, to 
examine the load-displacement behavior of the shell to 
determine the cause of convergence failure. A detailed 
discussion of the theory and implementation of SATANS can 
be found on pages 8 through 29 and in Appendices B, C and D 
Gr Ref. 3, and on pages. 8 through 28, pages 56 through 61 


and page 100 of Ref. 4. 


€. USER PREPARED SUBROUTINES 

Five small subroutines describing the shell geometry, 
the stiffnesses, the loading (physical and/or thermal), 
and the initial conditions (in the dynamic case) must be 


prepared by the user. Complete descriptions for the 
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onl of these subroutines can be found on pages 34 


gh 38 of Ref. 4. 
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LET. CONSTRUCTION OF SATANS=I 


Stein /Ref. 10/7 observes that a basic impediment to the 
widespread use of existing computer programs is the difficult 
and laborious process of assimilation of a program's basic 
theory, its implementation and its use by a prospective 
‘researcher. He also warns the ‘unwary analyst' of the dan- 
gers inherent in computer output, which is a ‘blizzard of 
numbers.' As a consequence, SATANS is continually being 
modified in an attempt to make it more usable. The most 
recent modifications made for this purpose are described in 


enesesection., This new version of SATANS is: called SATANS-I1. 


A. PACKAGING 

SATANS has been reorganized to allow 99% of the code to 
be prepackaged for repeated use by researchers. Compiled 
versions of all subroutines and main controlling routines 
can now be placed on a disk, in a data cell or on a magnetic 
tape, allowing the user to "submit" only the user written 
subroutines and inputidata. —A. complete listing of this 


storable package is given in Appendix E. 


B.° PLOTTING ROUTINE 

A great advantage is gained if the computer output can 
be presented in graphical format. Such plots, however, are 
often generated using mechanical or photo-electronic devices 


which are not universally usable. Consequently, SATANS has 
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been enhanced with an equipment-independent plot package 
utilizing self-contained subroutines and the high-speed 
line printer. 

The main program listed in Ref. 4 could not be modified 
EO, Incorporate the logic. and control for a graphical ability 
as it then exceeded the capability of the FORTRAN-H compil- 
er. Therefore, the main program of Ref. 4 has been replaced 
by a new main program and two subroutines described as fol- 
lows: 

SAGANS = Lhis is the new ‘“main.* it reads all data 
cards, prepares the title page of output and passes control 
to either of two basic controlling subroutines, 

SrATTC = Thas is the controlling subroutine for a static 
analysis. It conducts initialization, establishes geometry, 
stiffness, loading and boundary conditions, controls plot- 
ting of input information and controls the solution proce- 
dure. 

DYNAMC - This is the controlling subroutine for a dynam- 
ic analysis and it operates in a fashion similar to subrou- 
tane STATIC. 

The graphical ability is obtained by the introduction of 
myo controlling subroutines, PLOT1 and PLOT2, four plotting 
subroutines, PLOTIT, SCALIT, ROUND and DRAWIT, - and modifi- 
cation of subroutine OUTPUT for the selection of output data 


for plotting. A single data card with request information 


I the four plotting subroutines are slightly modified 
versions of those contained in Ref. 11 under the listing 
wo, ULPLOT, “ 
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is required to generate the plots. The plotting routines 
will handle a maximum of 100 stations. This routine is 


included in the listing given in Appendix E. 


C. MISCELLANEOUS ALTERATIONS 

The iteration lower limit has been raised to two itera- 
tions. This change was introduced because it was noted 
that immediately after a load step size reduction, several 
consecutive solutions were obtained in one iteration, a 
completely linear process. This was then followed by non- 
convergence and, in some instances, premature program termi- 
nation. Setting a minimum of two iterations automatically 
introduces nonlinearity into each load step and appears to 
solve this problem, as illustrated in Figure 4. 

A quadratic extrapolation technique was investigated 
using a static analysis of the highly nonlinear shallow 
spherical cap. No gain in solution accuracy was seen and 
a price was paid in program execution time due to approxi- 
mately 10% greater total iterations required, not to mention 
the increased storage space requirement. This phenomenon is 
exactly as observed by Stricklin et al /Ref. 127; hence the 
linear extrapolation technique has been retained. 

Since boundary loads are introduced through the boundary 
conditions and since the program considers the boundary 
conditions to be the same in every harmonic, changes indi- 
cated in Appendix B of Ref. 5 are permanently incorporated 


into the program. This alteration causes the applied 


25 




















NONCONVERGENCE AT THE NEXT 
LOAD STEP WOULD CAUSE PRE - 
MATURE PROGRAM TERMINATION 
HERE ( if the loadstep reduction 

limit has been reached) 


LAST SOLUTION PRIOR TO A LOAD 
\sTEP SIZE REDUCTION 


Loap 





ACTUAL SOLUTION PATH 
| ITERATION PATH 
— ———2 ITERATION MINIMUM PATH 


NOTES: 


|. HORIZONTAL SCALE EXPANDED 
2. NUMBERS INDICATE ITERATIONS 
REQUIRED FOR CONVERGENCE 











DEFLECTION 
PREMATURE PROGRAM TERMINATION 
FIGURE 4 


26 





- 5 4 
= a ae SS —— 








loading at a boundary to be taken in the zeroeth harmonic 
oniy, i.e., when n=0. 

All '"COMMON' statements have been ordered and separated 
and some slightly modified. This alteration allows easier 
reading and correlation of FORTRAN variables with specific 
subroutines. 

The output ability has been expanded to allow printed 
data (and plots) at up to 36 vice six meridians. 

riche. case of, a. Shell. with, an. initial..(or a final) pole, 
determination of the shear force at the pole in the n=1 har- 
monic has been deemed an inconsequential and lengthy computa- 
tion and has been removed. Appropriate changes have been 


made in subroutines INLPOL, FNLPOL and OUTPUT. 


Bee Usk OF SATANS-I 

instructions for preparation of input data. cards for 
SATANS-I are contained in Appendix C and the initiating main 
program and user written subroutines appear in Appendix D. 
Appendices D and E contain a complete listing of the 
computer program, SATANS-I. 

SATANS-I requires 206 k-bytes of main core when run on 
an) IBM-360/67 computer under OS/MVT. and compiled by the 
FORTRAN-H compiler at optimization level 2. Execution times 
are the same as for SATANS unless plotting is requested, 
and then is increased approximately 0.01 seconds for each 


100 point plot generated. 
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IV. CONSTRUCTION OF SATANS-II 


As discussed previously, a serious impediment to the 
consideration of initial imperfections by Dantone /Ref. 57 
was the requirement that the imperfections be symmetric 
about a datum meridian. This restriction results from the 
fact that all dependent variables were expanded in a single 
Sine or cosine series. The analysis presented in this 
section removes this restriction by introducing the use of 
full trigonometric series expansions for the dependent 
variables, and then modifies the resulting new relations 
where required for the introduction of initial imperfections. 
The computer program that incorporates these changes is a 


modification to SATANS-I, called SATANS-II. 


A. INTRODUCTION OF FULL TRIGONOMETRIC SERIES EXPANSIONS 
The original assumption by Ball /Ref. 37 was that a 
normal pressure load could be expressed in a Fourier cosine 


series as 


co 
q = (o.ne/e.) L p'™ cos né Ce) 
Meee OG = the weference jstress level, 
ho (= the reference thickness, 
E, = the reference Young's modulus, 
nm = jehe Pourver andex, or. harmonic, 
and p = the nondimensional Fourier coefficient of the 


normal pressure load. The normal displacement W 
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corresponding to this pressure was taken as 


W = (ac, /E,) ) w'™oos ne (2) 
0 


where a is the shell characteristic length and w is the 
nondimensional series coefficient of the normal displace- 
ment. For a completely arbitrary case, representing q and 


W with the complete Fourier series gives 


q= (oh. /E,) Ee p! (Moos no + ) p? (sin nd} (3) 
O 1 


and 


W = (ac./E,) ( } w! ‘Moos ne + ) w2'™) sin nd} (4) 
1 


1 indicates a previously existing coef- 


where superscript 
ficient and superscript * indicates a new coefficient. 
“Similar series expansions can be written for the other 
pressures, displacements, etc. 


Phe nonlinear terms associated with the full series are 


fidustrated by 62, 
s 


e2 = (0,/2,)*[ ( 1 42'cos no}*+ ( J 42) sin ne)? 
1 


ees ao 91 cos no} ( } ae né } ] (5) 
0 1 J 


where the %. are the nondimensional series coefficients of 
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the meridional rotation, given as 


guim) = “(-wi ))- 4 ea 
(6) 
jo os eve le r daa) 


Originally, in Refs. 3 and 4, the nonlinear term os was 


represented as 


5 + gi (n) y a(n 
o2 = [oe7EQ) 4 ) qe Cosrng) = 0 2/E, ) Ba) Cos mo (7) 


since 

Cos(n8)Cos(m8) = %Cos(n-m)68 + 4Cos(n+m) 6 
The pon term contains the products of all. the o's whose 
harmonic sum or difference equals n. There are now two 


additional terms in equation (5). These can be represented 


as 


(o27E,) > { d G2kO) san HO] =o. /E.. d 2M) cos né (8) 


and 


ZiOn7Eal (:) ¢1(™) cos ne} (1) o2() sin no} 
0 1 
(9) 
SY t/a. ) a UTE né 
1 


12 


where superscript indicates the product of a previously 


existing coefficient and a new coefficient. 
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Thus, 


ae aoe | } Bi ("cos ne + g2(M) cos ne 


(10) 


+} 822) sin no) 
1 


Similarly, the nonlinear term 0 (and also 6? since it is of 


the same form) is given by 


co 


93 = (90/Bo)*| (J 94 (sin ne)*+ ( J 02 cos na)? 


(11) 
+ 2( } 2 (ml sin no}( } Be ces noe } 
1 0 
which can be represented as 
o. 200 /ER uf d p1'n) cos no + ) p2'") cos née 
(12) 


foe) 


+ } pee sin no} 
1 


The nonlinear term o 2 is given by 


S) 
ae = [en /Ee ( ) $2 eos no.) (-) 64M) sin nd} 


1 


#() i lore mney fe) Cox nd} 
0 


o 


(1:3) 


ea 


oN sin ne} (-) 249) Sin nd} 
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~ri8g 


+ ( 5 


-C-1 
oO 
i) 
>| 
dp) 
}- 
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yo | 
@D 
—— 
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#2 cos no) 
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which can be represented as 


; 5 (n) 
$5% = F/E, (J ByaP) Sinyn® + } Be, 7 Cos nO 
(14) 


ze g2 (NM) o; 8 
) ae in née) 


Following this scheme, similar expressions can be obtained 
for the remaining nonlinear terms 2 Noe ®5No, PEN oy ®—5Ng,, 
®N, and Ng. 

All of the nonlinear terms fall into one of the three 
categories Cos(k@)Cos(£8@), Sin(k8)Sin(28) or Cos(k@)Sin(286). 
Since these three forms were already accounted for in 
SATANS, there was no need to develop any new routines for 
the full series expansions. It was only necessary to intro- 
duce each new nonlinear term into the appropriate existing 


summation in subroutines PHIBET and TEAETA. 


Be.) ENTRODUCTION OF INITIAL IMPERFECTIONS 
The initial imperfections only influence the nonlinear 
terms in the equations. The strain-displacement relation- 


ships used in this analysis are 


) 
M 


Ge W/ Re Ale + $7) 


s 
Ey =V /r + r“*U/r + W/R, 4 (06 = 67) (15) 
eo. vee Wee V7 oe e207) 
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where 


Oe aE UPR. 
o> =<W/ mt V/Rg (16) 
6 = &(v- + r°v/r - u'/r) 


and where ay, Eor € represent the reference surface 


s8 


strains, ® ) @ represent the reference surface rotations, 


SF age 
and superscript primes and dots represent partial differen- 
tiation with respect to s and 6, respectively. 

Consider a shell whose geometry differs slightly from a 
shell of revolution. For this imperfect shell, a deviation 
from a specific shape occurs in the form of a transverse 


displacement, W. This W produces no strain. Defining W* 


as a total displacement, 


Wwe =wt w (27) 


.and replacing W in equations (15) with W* of equation (17) 


gives 
e. + €. = U* + W*/R_ + 4 (0% + 0°) 
Ss Ss S 2S 
a f 2 2 
a ° ra Tk 1S o 
ep Gaya ea ruse WER at aos oe) (18) 
<= a 1 2 e a 2 
eo ee aca BVO ae eU fe = ero or on) 
where 
= uy 
ie W + U/R, 
(79) 


Wi re 
of W* /r + V/Re 
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and oy Eo and ae: are the pseudo strains due to W only. 
V ‘Ge 


When U = then We = W, “and ©] =e, “= Esp = 0. 
Substituting these relationships into equation (18), solving 


for € Eg, Esp and placing these terms in equations (16) 


Ss! 
lead to 
4 * ye q W L (Ww Z 
Eee O° + W*/RL + (o* +o} - W/R, - %(W") 
(20) 
2 i ae 
Bee VO /e + r*U/r + We7R, + 4(o% + o°} - W/Ry - (W") /r 
Ego TAO yeu ae Veo onog - Ww Jz) 


Expansion of equations (20) making use of equations (16), 


fe7) and (19) lead to 


= 4 I 2 2 ier =e 

ge WR (ol oO") = {Wo} 

cope /E ot EOU/E +. W/R, + 4(05) + 07). = {W"o,/r} (20) 
= ae 2 e = 4 f= Ty]? 77°e 

ome ot Ue VE + oO. = {WO + a Jey) 


where the terms enclosed in braces are the newly generated 


terms due to the imperfections, and 


tf =-o = {w= } 
s s 
(22) 
Che oe  larepe 
ey {w /r} 


The presence of imperfections has no effect on the changes 
of curvatures since the curvature-displacement relations are 


linear. 
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Expressing the imperfection partial derivatives -W~ and 
-W'/r as imperfection rotations ¢, and Py, the strain- 


displacement equations (21) can be expressed as 


m 
| 


4 2 2 — 
poe Us at W/o (Oe + (2c eh: 5 4 e?) 


Eg = Vi/r + r“U/r + W/R, + cia + (2a +507 t)S RES) 


20 B(V + U0" Jee ig {@ 9; + 6 9 }} 


: sg 


where the terms in braces are the newly generated terms. 
Imperfection rotations also affect the nonlinear term, ON. 
(and the others listed after equation (14) ). For these 
terms, the reference surface rotations $, and 5 must be 
replaced by o and 25 of equations (22). 

To illustrate the procedure which is used in the follow- 
ing section to introduce the initial imperfection rotations, 


consider the single series expansion, equation (7), repeated: 
here, 


co 


‘ a co 
2 = (0,/E,)?( } ¢1'™ Cos no)? 2 o,/E, } g1'™cos ne (7) 


Modifying this equation to include the effect of initial 


imperfections indicated by equations (23) results in 
{ Cc 
. a 2: 
ieee 20.0.) — (0,/e,)*| ( 61" cos ne)? 


Cos ne) | (24) 
y) 


Bs = 1(n 
= e/ Eo ee Vebs né 
0 
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A redefinition of B to include the cross-product term is 
now required, along with redefinitions of the other £8 terms. 
and the n terms in the other nonlinear term expressions. A 
redefinition of these terms alone is sufficient for the 
introduction of initial imperfections into a SATANS analysis 
where only the single series expansions are used. This was 


the thrust of the thesis by Dantone. 


C. COMBINATION OF FULL TRIGONOMETRIC SERIES EXPANSIONS AND 
INITIAL IMPERFECTIONS 


Consider that for every meridional rotation ®,, where 


®5 = (So/Eo}*( J ¢1@cos ne + 92) sin ne) (25) 


there exists an initial meridional imperfection rotation 


@., defined as 


(") sin no) (26) 


ae = (9/E,) * { ) $1 '™ cos ‘no a ) d 
The nonlinear term 2, previously given by equation (5), now 
becomes (o2 + 26.0.) , according to equation (23), where o< 
is given by equation (5) and @,¢, is the product of equations 
(25) and (26), as illustrated for the single series expan- 
sion in equation (24). All the nonlinear terms of equations 
(15) have been converted to those in equation (23) for the 
full series expansions and the resulting relations are 
given in Appendix A. The definitions of the 8 terms (given 


previously in Ref. 3) have been modified to include the 
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effects of both full series expansions and initial imperfec- 
tion presence. 

The remaining nonlinear terms (those listed after equa- 
tion (14) ) also require modification but are of simpler 
form since only the first power of the rotation is involved. 
A redefinition of the n terms to include the effects of full 
series expansions and appropriate initial imperfection rota- 
tions has also been made. 

In this manner, the validity of the strain-displacement 
equations (equations (5) of Ref. 3 and equations (15) of 
this report) and the equilibrium equations (equations (8) 

Of Ret. 3) is retained. Introduction of full series 
expansions and initial imperfections in the nonlinear terms 
was accomplished by a redefinition of the 8 and n terms, 
which are the nondimensional coefficients for the nonlinear 
terms in the strain-displacement and equilibrium equations, 


respectively. 


D. SUBROUTINE MODIFICATIONS 

Introduction of the full series into the program was 
implemented by use of negative index numbers for the new 
Sine or cosine portions of the expansions (the terms with 
the superscript 7). With the requirement that the axisym- 
metric harmonic (n=0) be included, the program was expanded 
to handle a maximum of 25 modes, i.e., the zeroeth harmonic 
(n=0) plus twelve other harmonics (positive and negative). 
Each positive index and each negative index will be referred 


to as a mode, i.e., the harmonic n=5 has two modes, -5 
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and +5. A typical correspondence between harmonic, mode 


and Fourier index is given in Table I. Harmonic zero (n=0) 


and any combination of up to twelve harmonics may be 


included in the 
Table I are the 
negative values 
the full series 
the first mode, 


fPeEest mode of a 


analysis. The positive values of NN in 
modes analyzed in SATANS and SATANS-I. The 
of NN are those modes added to incorporate 
expansions in SATANS-II. Note that n=0 is 
and that the negative value is always the 


pair of modes. This is a requirement. 


TABLE I 


HARMONIC, MODE AND FOURIER INDEX CORRESPONDENCE 











HARMONIC| MODE| FOURIER} FORTRAN ASSIGNMENT 
n INDEX STATEMENT IN SUBROUTINE 
PLOAD OR TLOAD 
0 i 0 








NN (1) =0 


-1 NN (2) =-1 
NN (3) =+1 
NN (4) =-2 
NN (5) =+2 
NN (6) =-5 
NN (7) =+5 
NN (24) =-18 


NN (25) =+18 
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Seven Seeoneenas of SATANS-I required changes to 
incorporate the full trigonometric series expansions and 
the initial imperfections. 

Subroutine PHIBET calculates the ¢'s and carries out 
the summation procedures for calculating the 8 nonlinear 
terms. Since extensive modification was required, the 
subroutine was altered such that those sections handling 
the full trigonometric expansions and imperfections are 
entered only if required, i.e., if a full series analysis 
is necessary. 

Subroutine TEAETA calculates the inplane forces and 
the n nonlinear terms. It was also extensively modified, 
again, utilizing new sections only if required. 

The subroutine MODES, a highly complex procedure for 
determining modal coupling required minimal alteration 
since each unique mode number is stored at every, or every 
other storage location, depending on use of SATANS-I or 
SATANS-II, respectively. Each unique harmonic is accessed, 
then, by an appropriate step size change in the controlling 
loops. 

The subroutine OUTPUT required the addition of a section 
to include the new half of the trigonometric series in the 
generation of a summed solution at a specified meridian. 

The subroutines POLE, INLPOL and FNLPOL compute the 
Eeiatron at an initial or a final pole. Due to the compu- 
tation of pole conditions based on n, which now has positive 


and negative values, these routines were modified 
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accordingly. Effectively, the only change in the equations 
presented in Appendix D of Ref. 3 involve a replacement of 
n by -n where appropriate. 

A complete listing of these seven modified subroutines 


is given in Appendix F. 


Be | SATANS-I1 

To modify SATANS-I to handle completely arbitrary 
loading and initial imperfections, the following steps 
must be performed. 

1. Replace subroutines PHIBET, TEAETA, MODES, OUTPUT, 
POLE, INLPOL and FNLPOL with those in Appendix F. 

2. Replace the affected COMMON statements in Appendices 
D and E with those listed in Appendix F. 

3. Add subroutine IMPERF, also listed in Appendix F. 
This subroutine contains the values of the imperfection 
rotation coefficients, $. and oy. Instructions for actual 
preparation of these coefficients will be discussed in de- 
tail in section V. If imperfections are not to be consider- 
ed, the o, and %% are never accessed. 

4. In Appendix E, make the following additions after 
the card indicated: 


Subroutine After Card Number Insert 


DYNAMC 00002270 CALL IMPERF 

STATIC 00002230 CALL IMPERF 

PMATRX 00001100 IF(N(MN) .LT.0)P(2,1,1J)=-1.0 
PMATRX 00002000 TF(N(MN) .LT.0)ZFP1(2,2)=1.0 
SATANS 00000380 READ (5,108) JUMP ,MPERFS 
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SATANS 00001400 108 FORMAT (215) 
5. In Appendix E, modify the following cards: 


Subroutine Card Number Should Read 


DYNAMC 00000800 DIMENSION VBAR(25) ,AVB(25) 
PMATRX 00000900 NN=IABS (N (MN) ) 
PMATRX 00001640 NN=IABS (N (MN) ) 


6. In Appendix C, modify the data preparation instruc- 
tions by adding a data card after card #4 according to the 
following format: 


Card Column Format Item Interpretation 


4.1 1-5 L5 JUMP 1 = single series expansion 
2 = double series expansion 
4.1 6-10 LS MPERFS 0 = no imperfections 
1 = imperfections 


An initial imperfection solution requires the use of full 
series in the analysis. Specification of the parameter 
JUMP indicates whether the solution will be as in SATANS-I 
or will utilize the full trigonometric series expansions of 
SATANS-II. Specification of the parameter MPERFS indicates 
whether the analysis will include initial imperfections. 
ieee that if MPERFS equals 1, JUMP must equal 2, but if 
JUMP equals 2, MPERFS may be either 1 or 0. 

we AGG imperfection data just prior to the last card of 
the data deck. This data is to be placed in arrays PHIXB 


ana, PHITB. 
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F. USE OF SATANS-II 

After completing the changes in section IV.E, SATANS-II 
exists aS a program capable of performing all the functions 
of SATANS-I with the following restrictions: 

1. The number of meridional stations cannot exceed 
Poo, (KMAX < 101). 

2. The number of modes cannot exceed 25, (MNMAX < 26). 

3. The product of meridional stations and modes cannot 
exceed 951, (KMAX*MNMAX < 952). 

4. Plotting ability is not available for the imperfec- 
iron input data. 

Specification of harmonics to be considered (accomplish- 
ed in subroutine PLOAD or TLOAD) must include the axisym- 
metric harmonic (n=0) and both negative and positive entries 
BOm any other harmonic, as indicated in Table I. 

The instructions for the preparation of the input data 
are given in Appendix C. SATANS-II requires 426 k-bytes of 
Main core on an IBM-360/67 computer when run under OS/MVT 
and when compiled by the FORTRAN-H compiler at optimization 
level 2. Execution times vary widely from four minutes for 
aeoo Station and 3 mode analysis, up to 30 -— 60 minutes for 


a 35 station and 25 mode analysis. 


G2. VALIDITY OF SATANS-IIL 

SALANS-LE was checked £06 1ts' ability to obtain a cor- 
rect solution using the full series expansions. This check 
was made uSing a clamped cylindrical shell subjected to 


pressure loadings as described in Figure 5. A solution was 
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C. ASYMMETRIC LOAD 


-- = SOLUTION COMPARISON POINTS 
SATANS -IL CHECK LOADINGS 
FIGURE: 
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obtained for the loading shown in Figure 5A with SATANS-I 
and the summed solution was printed at meridian 6 = 0°. The 
cylinder was then analyzed by SATANS-II for the following 
three orientations: 

1. Loading defined as in Figure 5A -—- a symmetric load. 

2. Loading defined as in Figure 5B - an anti-symmetric 
load. In this case, the summed solution was printed at 
Meridian 6 = 90°. 

3. Loading defined as in Figure 5C - a general load. 
In this case, the summed solution was printed at meridian 
Bee 45°. 

Load definitions for the four check runs appear in 


Table II. In each case, the solution obtained by SATANS-II 


TABLE II 


FOURIER COEFFICIENTS FOR CHECK LOADINGS 









FOURIER COEFFICIENTS OF APPLIED LOADING 
MODE | ------------------------------------------- 
SYMMETRIC | ANTI-SYMMETRIC | NONSYMMETRIC 









0.0 

















-1 0.0 -45.015 
+1 -63.66 -45.015 
-3 050 +15.005 
+3 +21.22 -15.005 
-5 0.0 +9.003 
+5 -12.73 +9.003 
-7 0.0 -6.431 







+9095 


Note: the coefficients are based on 1a) pS. 
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agreed with the SATANS-I solution” through four significant 


dagits. 


2This reference to solution includes the load-deflection 
history, the ultimate nonconvergence level, and the displace- 
ment, Lorce and moment distributions. 
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V. IMPERFECT SHELL ANALYSIS 


A. PREPARATION OF THE IMPERFECTION DATA 

With the developments of section IV implemented, 
SATANS-II requires only the Fourier coefficients of the 
imperfection data to conduct a nonlinear analysis. Recall 
from equation (26) that these Fourier coefficients, $. and 
% are determined by knowledge of the actual imperfection 
rotations, oF and @.. Recall also from equation (22) that 


6 


the imperfection rotations are given as 


¢ = -W- 
(27) 


oa. 
i 


-W'/r 


This means that the imperfection rotation Fourier coeffi- 
cients can be determined by knowledge of the initial normal 
deflection deviation, W. This imperfection information can 
be obtained from a mathematical description, a graphical 
plot or from actual experimental measurement. A simplified 
data generation flow chart is presented in Figure 6. 

1. Data Generation by Mathematical Description 

In a note on imperfect shallow spherical cap buck- 

ling, Kao /Ref. 13/7 analyzes a cap subjected to a mathemat- 
ically defined dimple imperfection existing over one-quarter 


of the shell. The dimple is given by 
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SIMPLIFIED IMPERFECTION DATA GENERATION 
SCHEME 


FIGURE:.6 
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W = 64(s/L}?(1-s/L}*hécos 26, -45°<@< 45° 


(28) 
Ww =O; 45°<6<325°9 


where L is the total length of the shell meridian and where 
s/L ranges from 0 at the pole to 1 at the edge, and 6 is an 
intensity factor, 0 <= 6 <= 1. The geometry of this problem 
is shown in Figure 7. Kao used a two-dimensional finite 
difference mesh of eight meridional and 16 circumferential 
stations in his analysis. A numerical model of, for example, 
35 meridional stations and 25 modes would yield a SATANS-II 
eatutcion of considerable accuracy. The discrete initial 
normal deflection deviations W at each mesh point need not 
be computed, as partial differentiations of equation (28) 
give expressions directly solvable for discrete values of 

a and ®, for each s and 0. Imperfection rotations at loca- 
tions exterior to the affected area (=45° <6 < 45°) are 
zero. To illustrate, the meridional imperfection rotations 
that exist along the line ABCDEA of Figure 7 would appear 

as curve A of Figure 8. The discrete values generating 

this curve can easily be expanded into a full Fourier series 


(m) inserted into “SATANS-I1 


and the resulting coefficients $. 
for the analysis. 

The imperfection of Figure 7 is symmetric about the 
meridian @ = 0°. As an example of a more complex imperfec- 
tion, another dimple of the same form could be located in 


the region defined by 100° < 6 < 190°. In this asymmetric 


case, the meridional imperfection rotations along line 


48 


‘> 
<<" 
~ 


aed f toes eld: Se Bh esos ae oF 


x 7 


ii es Rites 4 A4-se r “eye me pass 






















, ‘oh { Ab 
4 cee: set - is) Oo ee ih este “Pig 
en GRIT LE WY Een 1d -s Cr an ahh A 
j | ; ae Lee 7h a ot. Sioks ‘iG, iene dd. 
: a 
7 10 faba CAS tae ¥ BLUYL GAG By Ke, th re - ’ 
4 in ek ae = pe aH shiver 8 TTR nine. 2 e2h repose 
sul iY ihe VIAL TOIEe Wide "ge ae a tw 
(Lod deem, pVeg, 2H Ww pre toaiven. ee a (39% 
OL Bolded Ro Sota +70°7n 139m 4008389 ee , tet ge r 
ea (seams Apres =guy gre at a @thient " ky swat fe oe ’ Ts 
achtwice 2Oki Ge oem x 1 pte @, ae y a on ne 
alt : or vo) SR cE Tite? oe. OTS a Tad tT? xo 
TE? Palettes 2a Lane bi 165 gtd (oh ess agit: of =a 
/ tok’. te IR, que, ji ad’ enole Setaem 
hati ist eouts ihe 4 iain k Gy opt % osegtt io & eFaee 
i! rus ge’ be n i’ fai 14itNS ‘oe bees" fats SVieat 


An Bangin hal 3 Wed Sab Ra Meee ends lyataans sa 


eptaylonn ni 


tee al K anys 2B:  eaktoe® Taaid Ts a 


wT age) »~e hams =a ih te “ak mint! in aA e wy s 
| Soranhd, oe Bia “ar aoa eet Aiea ae 
: = 


Ahern Syne eons oN Parti a) ¥ ae. , 4 te 
; ' ee Re y _ i oat 
o% ett J [ NOt fer) oe Wh i rs 6 ae. ho ‘- 4 poet: ae ‘ oe - 








@=45° 
s) 
a ' 
=190° | ak 
8-190". ce 
O= 325° 
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IMPERFECT QUARTER WITH 
DIMPLE IMPERFECTION GIVEN BY 


W/h= 64x°(1-xp dCos20 
WHERE X= sf 


MATHEMATICAL IMPERFECTION GENERATION 
FIGURE 7 
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A.SIMPLE DIMPLE MERIDIONAL IMPERFECTICY ROTATION 





B. COMPLEX DIMPLE MERIDIONAL IMPERFECTION ROTATION 


MERIDIONAL IMPERFECTION ROTATION 
REPRESENTATION 


FiGURE (8 
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ABCDEA of Figure 7 would now appear as curve B of Figure 8. 
The discrete values generating this curve can also be easily 
expanded into the full trigonometric series and its coeffi- 
cients inserted into SATANS-II. 
2. Data Generation by Graphical Plots 

Graphical representations of W can be easily con- 
structed and the patterns are near limitless in design or 
location. Theoretical or actual creases, furrows or dents 
of any shape need only be graphed, the discrete imperfection 
displacements determined, the discrete imperfection rota- 
tions ips and on computed via a finite differencing tech- 
nique, the Fourier expansions performed and the resulting 
eoetficients inserted into SATANS-II. 

3. Data Generation by Experimental Measurement 

Arbocz and Babcock /Ref. 6/ have developed a means 
eeemapping the topography of circular cylindrical shells, 
and can thus describe the deviations of an actual shell. 
The procedure, described in detail in Ref. 6, produces 
imperfection data in the form of nondimensional normal 
displacements from a perfectly circular cylinder. For the 
sample shells considered here, 49 data points were taken 
around each circumference and 29 equally-spaced circumfer- 
ential traces were made. The first and last trace were 
0.64 inches from the ends of the shell and the spacing 
between traces was 0.24 inches. Pertinent geometrical and 
material data concerning the shells of interest here is 


Contained in Table III1. 
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use by SATANS-II, 


dimensionalized by the wall thickness. 


TABLE III 


IMPERFECT SHELL OVERALL CHARACTERISTICS 








SHELL DESIGNATION 












Radius 4.00" 4.00" 


Length 8.00" 8. 00" 

“Young's Modulus x 10°, & [15.10 psi [15.80 psi 
‘Wall thickness, h —*‘|- 0.00444" | 0.00437" 
“Number of axial traces | 29 4| 29 
‘Axial increment, A || 0.24" ~— | 0.24" 


ee ee me es me ee ee ee en ee ee ee ei oe me ee eee ee ee ee es oe ee ee ee ee ee ee ee ee 


Number of data points 


To prepare the imperfection data on shell A-14 for 


To retain the 


original axial increment of 0.24 inches between circumfer- 


the nondimensional deviations were first 


ential traces and to make the data compatible with a 35- 


station meridian, the shell length was increased by 0.08 


inches at each end, or a total of 2% 


(See Figure 9). 


required three additional traces at each end of the shell 


corresponding to stations l, 


2 ange 3 


and, 33:, 


S34. .and Si5% 


Bata for “the additional six circumferential traces were 


generated via the modified linear extrapolation procedure 


mblbustrated in Figure 10. 


the meridional and circumferential rotations were computed 


Reveach othe, i p7Ssdata points, 
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(REF. 6) (SATANS- IZ) 


NOTE: Extrapolation required for data 
at primed stations 


CIRCUMFERENTIAL TRACE NUMBER AND 
MERIDIONAL STATION NUMBER CORRESPONDENCE 


FIGURE 9 
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@ = KNOWN DATA POINT (REF.6) 

H=TWO STATION AVERAGE 

4 =EXTRAPOLATED DATA POINT 

Wiz IMPERFECTION DISPLACEMENT 
OF THE ith STATION AND THE 
jth MERIDIAN 


NOTE: NUMBERS CORRESPOND TO 
SATANS-IL MERIDIONAL STATIONS 


IMPERFECTION: FREE MERIDIAN 
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using central finite differencing in the appropriate direc- 
tion. After nondimensionalization according to the 
SATANS-II scheme, each circumferential vector of 48 
meridional and 48 circumferential imperfection rotations 

was expanded via a periodogram into a Fourier series accord- 


zag to 


Se = 70) , #*. <c(n) ;’ <s(n) 
o/ogt, - o. oF d Tee eOse ne) Gan \ sumene,, me29) 
1 


and 


24 
(E,/o,)%, = Go el Gon Gos le Ge cin nd (39) 
1 


0 1 
Aemthis point, the Fourier coefficients 9. and $4 are -of 
preeper form for direct insertion into SATANS-II. The 
computer program which performed the actual imperfection 
data processing for shell A-14 is given in Appendix G. 
4. Selection of Important Harmonics 
Since the number of meridional stations was 35, the 


program maximum of 25 modes can be used in the analysis.° 
Note, however, that in equations (29) and (30), there are 
49 modes describing the imperfections. Therefore, a selec- 
tion of a maximum of 25 imperfection modes must be made. 
The harmonics to be considered in the analysis should in- 
eeuder those critical to the shell and those significant 


in the imperfection. The description of the harmonic 


3Recall the restrictions given in section IV.F. 
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selection process for shell A-14 is given in a later 


section. 


B. BOUNDARY CONDITIONS 

The testing configuration used by Arbocz and Babcock 
on the axially loaded cylinder is shown in Figure ll. The 
ends of the shell are attached to an end ring and a load 
cell with cerrelow, a low melting-point alloy. The ends of 


the shell are believed to be clamped so that 


iS 8 = Vi ="0 (31) 


The axial load on the test specimen was applied by means of 
four screws located at the corners of the square end plates. 
At each load increment, the screws were adjusted to give as 
uniform a circumferential load distribution as possible. 

In all versions of SATANS, the boundary conditions are 


represented by the matrix equation 


N U 
Na * V 
ff +/A7 = {2} (32) 
Q W 
Ss 
 . M. 


Boundary conditions can be specified on N, or U, N CEI, 


S68 
QO, or W, and %, or M, by prescribing input values for the 
matrices / 2 / and / A 7. There is one equation (32) for 


each boundary of the shell. In their present form, SATANS-I 
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and SATANS-II assume the column matrix {2%} as zero for all 
harmonics except the axisymmetric harmonic. 

As was mentioned earlier, boundary condition discrepancy 
between mathematical model and physical experiment is cited 
aS a major contributor to thin shell buckling load disagree- 
ment. For this reason, three interpretations of Figure ll 
were made. The cylinder was represented as experiencing an 
axisymmetrically applied force N, on both ends (condition A), 
an axisymmetrically applied force N. on one end and a 
meridional deflection restraint on U at the other (condi- 
tion B), and an axisymmetric incremental end shortening U 
at one end with meridional deflection U restrained at the 
other (condition Owe These three conditions were investi- 
gated on the clamped cylinder, and on a simply supported 
cylinder for further comparisons and validity checks. 
Mathematical assignments for boundary conditions A, B= and 


Gane shown in Table IV. 


C. ANALYSIS 
me. Duckling of the Thin’ Perfectly=-Circular ‘Cylinder 
A necessary part of this attempt to develop a pro- 
GeamerOr predicting the buckling load of imperfect shells is 
the establishment of the program's ability to predict the 
buckling load of perfect shells. Consequently, this section 


compares the results of a SATANS-II analysis of the thin 


foe is believed that boundary condition C best repre- 
Sencs, the actual boundary conditions in Figure ll. 
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cylinder without imperfections with the bifurcation buckling 
analysis presented by Almroth /Ref. 147. Almroth's study 
deals with the influence of edge constraints on the critical 
oad of cylindrical shells of the type used by Arbocz and 
Babcock. Equation (33) is given by Almroth for the axisym- 
metric displacement shape of the simply supported cylinder 


(Boundary condition A) prior to buckling. 


W = x (2vyN. + p) (2 ++A,Sin(a,s)Sinh(a,s) 
(33) 
+ A,Cos(a,s)Cosh(a,s) } i 


where 


Clon K(1 = N.)? 
K =(2y) 7? 


/as/ (2 = 2) \2 


2 
| 


Ny =-N,/(2yEh) 


(pr) / (Eh) 


O | 
M 
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({1 - N2}¥%sin(a,)Sinh(a,) + N.Cos(a,)Cosh(a,) ) 


A, = > 
(1 - N2) %(Cosh? (a,) + Sin?’ (a,)} 
({1 - N2}%cos(a,)Cosh(a,) - NSin(a,)Sinh(a,)) 
A2 —) 
(1 - N<) ?(Cosh* (a2) - Sin*(a.)} 
@, = (a,L)/(2r) 


G2 = (a,L)/(2r) 


and where No is the axial load, p is a uniform normal pres- 
sure, h is the thickness, E is Young's Modulus, v is 
Poisson's ratio, L is the length and r is the radius. This 
solution is referred to as the consistent prebuckling solu- 
tion, as opposed to the membrane or classical buckling solu- 
tion where W is a constant along the meridian. Equation 

' (33) has been solved for shell A-14 and W is plotted as a 
function of s in Figure 12 for an Ne = 0.21904. Also shown 
in Figure 12 is a plot of W versus s obtained from SATANS-II 
for the same N. for both 35 and 69 meridional station anal- 
yses. Note in Figure 12 that even though only a few sta- 
tions lie in the boundary zone, the trend is towards the 
correct solution. The displacement of the 35 station solu- 
tion, however, is only a rough approximation of the actual 
state. 


Almroth has determined the minimum bifurcation 


buckling load, Nsoy, and the corresponding value for n for a 
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range of parameters L/r and r/h, using both the consistent 
prebuckling displacement W given by equation (33) and the 
membrane or classical prebuckling solution, where W is 
constant, in the Donnell equations. For the simply support- 
ed cylinder with an r/h of the same order of magnitude as 
Sao A-14.and an L/r = 0.5, the critical. lead is Ng p= 0-863 


(N =39.5 lb/in) with n=27 for the consistent analysis, and 


Scr 


N Or ttn 


i ae lb/in) with n=27 for the membrane pre- 


Sc 
buckling analysis. Thus, the consideration of the axisym- 
metric edge effect, which is analogous to an imperfection, 
reduces the buckling load by 14.1%. Parametric studies 
conducted by Almroth show usage of from 100 to 800 merid- 
ional stations over the half-length of the shell. 

As shown in Ref. 5, use of the membrane prebuckling 
displacement in Sanders' equations for the simply supported 


cylinder simplifies the problem to an eigenvalue problem. 


The buckling mode is 


W = 6 Sin(mtTs/L) Cos (né) (34) 


The equations for the matrix elements in the eigenvalue 
problem are given in Appendix B. A plot of the critical 
Values of m and n is given in Figure 13. The minimum buck- 
ling load occurs when n=27 and m=15, and corresponds to an 
Nscr Of 475 lb/in. 

SATANS' ability to predict the minimum bifurcation 
load of a shell for a specified value of n has been demon- 


Strated by Stilwell /Ref. 15/7. SATANS-II analyses of shell 
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CRITICAL BUCKLING MODES FROM EIGENVALUE 
ANALYSIS, SHELL A-14 


FIGURE 13 
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A-14 (simple support, boundary condition A) were conducted 
over a wide range of n for both the consistent and membrane 
prebuckle cases. Figures 14 and 15 give plots of the crit- 
ical buckling load versus the number of meridional stations 
for various values of n. The dashed parts of the curves 

are estimated extrapolations. It is clear that as the 
number of stations is increased, SATANS' critical buckling 
load decreases and the associated harmonic increases. While 
the n=27 curve predicts the highest buckling load over the 
range of stations used in both cases, extrapolation of 

these curves indicates that the n=27 curve would predict 

the lowest buckling load if the number of stations exceeded 
200 to 300. Note that the membrane prebuckle critical 
buckling curves of Figure 15 seem to flatten much more 
quickly than for the consistent case. This appears to 
Gontarm the occurence of a reduction in critical buckling 
load due to the presence of an axisymmetric edge restraint, 
as noted by Almroth. Figures 16 and 17, which present plots 
of the critical buckling load versus critical buckling mode 
Gemowel? A-14 for the 35, 69 and 100 station analyses, also 
predict this lowering of buckling load due to edge restraint. 
the 'fanning' of curves in Figure 16 indicate an increase 

in the number of meridional stations will produce a lesser 
buckling load decrease than in Figure 17. These figures 
also predict a lower critical buckling load and a higher 


critical buckling mode with increasing meridional stations. 
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Figures 14 through 17 indicate that, given enough 

peations, SATANS-II can predict a critical buckling load 
and mode for the thin-walled shell that agrees with 
Almroth's prediction. However, for the 35-station analysis, 
which is the analysis to be used for the imperfect shells, 
mieeecritical load of 51.1 Ib/in and harmonic of 16 are 
Significantly different from the correct values. Asa 
result of the coarse mesh used by SATANS-II in the following 
analysis, the critical buckling loads for the imperfect 
shells will be higher than the experimental values. 
Nevertheless, it will be shown that the introduction of a 
totally arbitrary, naturally-occurring imperfection pattern 
on the imperfection-free mathematical model of shell A-14 
Besmlts in a Significant reduction in critical buckling 
load. This is not only sufficient to demonstrate the 
capability of the program to perform the analyses, but it 
also provides an estimate of the degree of imperfection 
- sensitivity of the shell. It will be shown that application 
of this estimate of sensitivity to Almroth's buckling load 
provides an accurate prediction of the buckling load of the 
imperfect shell. 
7A ieeodue lion of Imperfections 

When prepared as indicated in section V.A.2, the 
imperfection rotation Fourier coefficients are in a form for 
direct introduction into SATANS-II. Figure 18 shows the 
applicability of the periodogram method of expansion of the 


imperfection data according to equations (29) and (30). 


70 





‘1¢} 










cork, 13. OSS] sie : APES 


iste. bel tex-nins Bea 35% Sham) 


45), pate = ‘tise 


f 125 te) S07. s+ 
= 
yond o}3fivh et ~Asee 
ages = . : 
PY aA = 
- noe 3 ; 4 ol 1 ry" et ‘oe. | a 
“ 3 
{ omer 2G 
I é S| de a i eo & 7] 
eae a) a ets 438s re ve 
| ‘ , 
snes anyt-notssetseqms 


vite ore Yind Joc ee ee 
= 4 - é q 
= : 
okisa ot maspord Sas ne Vp tw 
Sit. t 43 Sm. 2k A& 264.96 
to ‘met Pie at pied) ens ae ytivio 
. .4 vhiyittenoe Jo eiantses 
Fy ° , 


me 


sanztostteqgas: 3S noksniboztay) : 4 
a.¥ soljose 0 Badeatial ee ca agian aa 


-* o¢@ esmetol steve a. anvaeeees ise 


2503 SiVid -e sl2adnois 
=) 
@ SA d - SNOILVLNAS3Yd3y NOILVLOY NOILOSSAYSdWI IWILNSAYSSWNOYIO 





NOILVINAS 3¥d3ay¥ YSINNOS DJINOWYVH G3L935135S €| 
NOILVLINSSSAYd3SY Y3IYNOA DINOWYVH Ge 


VeLV GC) Sea © 


6091 3-08 oO OlNORIVLS) JeSOVELlevl=Vv TWaHs 
NOILVYNDIANOD NOILDSSYSdW! 31dWVS 


c 


eusoq xX suDIpDs 
NOILVLOY TIVWILNSYSSWNIYIS 





ab 


=> - © 


4 


+ 
Y 


* at 


i 


r. 
¢ { 


sS9R Seg 





é 
% 
a at 
: 
ae 
ad > | 








sia 


i Z es - = f= - = 
he - a = 
i] > a 


id 


240 23LV1d gi aUYNDIS 
@ sA°O -NOILVLNSSSUd3Y NOILVLOY NOILOSSUSdW! IWILNSYSSWNOYIO 


NOILVLNSS3Y¥d3Y YSINNOS DJINOWYVH G3L9373S FI 
NOILVLN3S3Yd3Y YSIYNOS DINOWNVH SZ 
Vivd 9'43uO 
o09E 705.081 ‘O! NOILVWLS ‘2 3OvyL ‘bi-vy IWSHS 
NOILVYNOISNOD NOILDSZSNSdWI AIdGWVS 








f 
O 
.9 
ey/s/e°Q xX suUDIpDS b 


NOPEVICd WILNSYSSNNDYID 


72 


Tn! 


“hb 


\ © 
a — 
24 eae Pd 
4 
“~ 
4 
= 1a 
25 
oe” 
ay i 
=e 
- 
' T 
> oy 
; tm 
7 S| 
oy ee u 
ty ok 
= 7 
a ve ' 
> te > 
ma a 
® 
i] 
i] 


Ftv 


LOBE Hi 


S9Rere Wweeee Ft lion 





Curve A is the graph of a solution of equation (30) pre- 
senting the circumferential imperfection rotation at 
station 10 plotted against @. It is seen to pass smoothly 
through all the data points. 

Only the zeroeth harmonic (n=0) and twelve others 
can be accommodated by SATANS-II. In the case where the 
shell response and the imperfection are limited to a few 
major harmonics, all important harmonics can be included in 
the analysis. For example, the A=6 shallow spherical cap 
has only the n=2 harmonic as the major harmonic of response, 
and the imperfection curves of Figure 8 require only a few 
Fourier harmonics for adequately accurate representation. 
Also, the few required harmonics for imperfection represen- 
eatvon are in the vicinity of the n=2 harmonic. This pro- 
blem is will within the limits of SATANS-II. 

The thin-walled cylinder, however, responds in a 
large number of harmonics. This can be seen in Figure 16 
-- many harmonics produce critical buckling loads close to 
the value of the minimum critical buckling load. Addition- 
ally, the imperfection requires many harmonics. Thus, as 
memtroned in section V.A.4, a difficult selection of harmon- 
ics for use in the analysis must be made. Six harmonics 
were selected on the basis of the shell response and six 
harmonics were selected on the basis of the imperfection 
Ghnaracteristics. Of course, all 13 selected harmonics will 
contain both the imperfection characteristics and the shell 


response characteristics. 
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For the shell response harmonic selection, the 
critical buckling loads of shell A-14 for specified values 
of n, discussed in the previous section, are presented in 
Table V. Values were assigned to each harmonic according 
to the legend of Table V and six were selected based on 
these assignments. 

For the imperfection harmonic selection, the sums of 
the absolute values of the Fourier coefficients for both 
sine and cosine terms of both meridional and circumferential 
imperfection rotations were tabulated and are presented in 
Table VI. Values were again assigned to each harmonic ac- 
cording to the legend of Table V, and six additional harmon- 
ics were selected. 

The complete harmonic selection process for shell 
A-14 was accomplished by a comparison of the values assigned 
to each harmonic for both the imperfection and the shell. 
Table VII presents this comparison along with the actual 
harmonics selected for the analysis. Note that no overlap 
of important harmonics occurred, as the imperfection re- 
quired harmonics seven and below and the shell required 
harmonics ten and above. The harmonics used with shell A-13 
were n=0 and n=2,3,4,6,8,9 for the imperfection and n=10,1l, 
12,13,14,15 for the shell response with a simple support, 
and n=13,14,15,16,17,18 for the shell response with a fixed 
support. 

Once’ the) Specific harmonics for inclusion. in the 


analysis are determined, the appropriate 6, and %, 
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TABLE V 


CRITICAL BUCKLING LOADS OF INDIVIDUAL HARMONICS 












FIXED SUPPORT SIMPLE SUPPORT 
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The trigger in each harmonic is an internal 
pressure of 107’ psi. 


Each analysis included the indicated harmonic 
plus n=0. 


Value assignments are according to the following 
legend: SYMBOL | VALUE 
— absolutely necessary 
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TABLE VI 


RELATIVE IMPERFECTION ROTATION COEFFICIENT SIZES, SHELL A-14 


* 
* 
+ 
=. 
+ 
* 
+ 
+ 
+ 





Notes: -1. Numbers represent the sums of the absolute values 
of the imperfection rotation coefficients over 
the entire shell. 


2. Value assignments were made according to the 
legend of Table V. 
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TABLE VII 


HARMONIC SELECTION, SHELL A-14 


SHELL SELECTED 
RESPONSE HARMONICS 





FIXED SIMPLE FIXED SIMPLE 
SUPPORT | SUPPORT SUPPORT SUPPORT 


X X X 
X X X 
Pe X X x 
op eae X X Xx 
+ | * Xx X X 
+ | * X X X 
~ - 
* | 4 Xx X X 
+} = 
+ + * - X 
-_ == * = X : 
+ + * + x x 
- | + * x * X X X 
- | - * * X Xx Xx 
SS ; - X X X 
- * *K* S% X 
= = * * x X 
= = * * x 
- + * 
- = 
- = 


Harmonic zero is required by SATANS-II. 
Value assignments are taken from Tables V and VI. 


'x' indicates a selected harmonic. 
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coefficients (obtained from the periodogram) must be loaded 
into the SATANS-II arrays PHIXB and PHITB (performed in 
subroutine IMPERF) in the following order: 

i. ‘The £irste KMAX’ locations Contain’ the 95's (or $,'s) 
of mode 0. 

2. The second KMAX locations contain the $,'S (or $g'Ss) 
Of mode ~—i. 

3¢ The third KMAX locations contain the $.'s (or $9's) 
of mode +i. 

a. ©The fourth KMAX locations contain the $,'s (Or $,'s) 
of mode -j. 

Doeehee cetera, unveil all modes are-entered. 
KMAX is the FORTRAN name for the number of meridional sta- 
tions to be included in the analysis. One special case must 
be considered where the shell has an initial or a final 
pole. In this case, rotations are zero in all but the n=l 
Heammonic. If this harmonic is included in the analysis, 
the imperfection rotation coefficients at the first (or 
last) station should all be set to zero except for the 
Seertficients in the n=l harmonic. In the case of an initial 
pole, the coefficient associated with the -l index should be 
placed in the first location of the second set of KMAX loca- 
tions in arrays PHIXB and PHITB, and the coefficient assoc- 
lated with the +l index should be placed in the first 
m@cation of the third set of KMAX locations. For a final 
pole, the coefficient associated with the -l index should be 


Beacedsdniithe Liest location of the next to last set of KMAX 
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locations, and the coefficient associated with the +l index 
puould theeplaced im the first location of ithe Last ‘set of 
KMAX locations. 

To illustrate the effect of the selection of har- 
“monics on the imperfections, equation (30) was solved using 
only the 13 selected harmonics. The resulting circumferen- 
tial imperfection rotation representation at station 10 
appears as curve B in Figure 18. To determine the effect of 
the inability to include all nine harmonics deemed important 
to the shell response with a fixed support, an alternate 
selection of harmonics was made for this case (see Table 
VII). A complete SATANS-II analysis was performed with 
this alternate selection and the results were found to 
differ from those using the primary harmonic selection by 
approximately 2%. 

BS. Results 

‘* Six analyses (boundary conditions A, B and C with 
clamped and simple supports) were made for each shell where 
the imperfections were reduced by a factor of 1000. These 
‘minute’ imperfections acted only as triggers for the modal 
deflections and the resulting solutions gave the shell's 
eensistent critical bifurcation buckling loads. These 
results are presented in Table VIII under the category of 
Minute imperfections. Selected cases were analyzed where 
the minute imperfection triggers were replaced by a 10° ’psi 


internal pressure trigger and the results were identical. 
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TABLE VIII 


BUCKLING LOADS, MINUTE VERSUS FULL IMPERFECTIONS 








SHELL A-14 SHELL A=-13 


SUPPORT| IMPERFEC-|BOUNDARY |BUCKLING| REDUC- | BUCKLING | REDUC- 
TION LOAD TEON 
(3) 


(1b/in) (3) 





@ 

L Minute B 

A @ 

M 

P 

E A - 

D Full B O55 
€ 22.10 
A 

S Minute B 

af Ec 

M 

P 

L A 26530 

E Full B 24.90 
(a 22.80 


Rersutfticient data. 


The same six analyses for each shell were also made 
using the full imperfections. Results for the final load at 
which a converged solution could be obtained are also shown 
in Table VIII along with the percent reduction in the buckling 
load for each case. 

Since the full imperfection is introduced at the 
first load step as a pseudo load, the program requires con- 


Siderably more iterations than normal to obtain the first 
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solution. The subsequent solution usually proceeds ina 
completely normal fashion. However, convergence difficul- 
ties may be encountered in cases where the imperfection 
level is high. This was the case with shell A-13, as some 
of its imperfection deviations exceeded 3% shel wall thick- 
nesses, and fluctuated over 6 thicknesses in a @ arc of only 
53°. For shell A-13, a 70% initial imperfection level was 
the maximum for which a converged solution could be obtained 
on the first load step. Consequently, buckling loads were 
computed with a maximum of 80%, 85%, 90% and 95% of the full 
imperfection data. For the 80% imperfection, 70% was intro- 
duced at load step one and 10% was added at load step two. 
For the three other solutions, 5% increments were added at 
various load steps. Since SATANS uses the total value of 
the variables and not the incremental values, each inter- 
mediate solution is correct for its given load and imper- 
fection level, but the load-deflection characteristics are 
invalid. An imperfection level of 95% was the absolute 
maximum for which a converged solution could be obtained for 
for shell A-13, and a short 5% extrapolation was made to ob- 
tain the data presented in Table VIII for this Enea ge This 
procedure is illustrated for the shell with a simple support 
in Figure 19. A comparison of load steps and iterations for 
both shells under these conditions is presented in Table IX. 
>For Shell A-13 with clamped boundary condition A, a con- 
verged solution was not obtained for an imperfection level 
Boeove 85%. This condition is therefore not included in the 


results. A similar problem was encountered by Dantone /Ref. 
57 with this boundary condition. 
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SIMPLE SUPPORT 
SHELL A-I3 


—= == == EXIBAPOLATIONS 10° JHE I0O% 
IMPERFECTION LEVEL 
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CRITICAL BUCKLING LOAD vs IMPERFECTION LEVEL 
FIGURE 19 
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TABLE IX 


LOAD STEPS VERSUS ITERATIONS 


SHELL A-14 SHELL A-13 






LOAD|INCREMENTAL | ITERATIONS |LOAD| INCREMENTAL | ITERATIONS 
STEP | IMPERFECTION | REQUIRED STEP | IMPERFECTION | REQUIRED 






























ADDITION ADDITION 
L 100% 7 L 703 24 
2 2 2 10% 6 
3 2 3 53 16 
4 3 4 12 
5 3 5 9 
6 Ss 6 3 
fi 3 q) 8 
8 3 8 6 
9 3 9 5% 19 
10 3 10 3 
ei 3 Les 53 60 
2 4 12 58 
13 4 i es 18 
14 4 14 16 
5 5 5 16 
16 7 16 12 
£7 10 
18 16 1S8t load step size reduction 
1St load step size reduction| 17 | 4 
19 3 2nd load step size reduction 
EXO) 2 
18 3 
294d load step size reduction] 19 2 
2a | 3 3°4 load step size reduction 
3rd load step size reduction termination 
22 | 2 (Poy=40.6 1b/in) 
termination 
(Por=46.0 lb/in) 
Hoce: This data corresponds to boundary condition B, simple 


support. 
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D. DISCUSSION 

The SATANS-II investigation of the buckling load of the 
perfect cylinder (discussed in section V.C.1) indicates that 
a numerical analysis using 35 stations results in a higher 
predicted buckling load than the correct analytical buckling 
load. Consequently, the SATANS-II results for the predic- 
tion of the critical buckling load of the imperfect cylinder 
are expected to be high. Table X compares the results of 
the SATANS-II analysis (column A) with the actual experi- 
mental results of Ref. 6 (column E) and shows this to be 
the case. The use of more stations to lower the predicted 
buckling load can be accomplished only by either decreasing 
the number of harmonics considered in the analysis or by 
increasing overall program size. Both of these choices are 
impractical at the present time since the number of harmon- 
ics already in the analysis appears to be low considering 
the number required for shell response and the imperfection 
data, and the present program core storage requirement is 
large. 

-It is possible, however, to apply the SATANS-II results 
for the percent reduction in the buckling load due to the 
imperfections (see Table VIII and Column B, Table X) to a 
more accurate buckling load of the perfect shell from an 
analysis like that of Almroth to get a more accurate estimate 
of the buckling load of the imperfect shell. When the per= 
cent reduction from the consistent minute imperfection 


buckling load to the consistent full imperfection buckling 
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load, determined by SATANS-II, is applied to the accurate 
consistent bifurcation buckling load of Almroth (column C, 
Gable X), the result is a prediction of the actual buckling 
load of the imperfect shell. These imperfect shell predict- 
ed buckling loads are shown in column D of Table X and are 
seen to be very close to the actual buckling loads of the 
imperfect shells, given in Column E of Table X. 

The percent reductions in buckling loads determined 
by SATANS-II are considered to be applicable to the accurate 
consistent bifurcation buckling loads of Almroth since the 
imperfection characteristics are seen to be concentrated in 
the lower ten harmonics (see Table VI), where shell response 
does not appear to be largely affected by the number of 


meridional stations in the analysis (see Figures 14 and 15). 


fee PROPOSED ANALYSIS PROCEDURE 

When the SATANS-II predicted buckling load of the perfect 
shell is higher than the analytical prediction, as is the 
case with the extremely thin-walled cylinder, the following 
procedure is proposed. 

1. Determine the consistent bifurcation buckling load 
as close as possible. 

2. Determine the buckling load with SATANS-II using 
Minute imperfections. 

3. Determine the buckling load with SATANS-II using 
full imperfections and compute the percent reduction from 
step 2. 


4.- Apply this reduction to the buckling load of step l. 
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This procedure is expected to work only when the result 


oO: Step 2 is reasonably close to that of step 1. 
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VI. CONCLUSIONS 


A. SATANS-I 

The computer program SATANS-I has been constructed from 
SATANS. It is a more user-oriented version of the program 
by Ball, whose primary modification is the introduction of 


an equipment-independent graphical output ability. 


B. SATANS-II 

SATANS-II is a further modification of SATANS where the 
mestriction that all applied loading and initial conditions 
be symmetrical about some meridional plane has been removed. 
In addition, SATANS-II possesses the capability of including 
experimentally-determined naturally-occurring imperfections 
or malformations induced during fabrication, mathematically 
described imperfections of any type, frequency or location, 
or graphically described imperfections due to damage suffer- 
ed after fabrication. The imperfections are introduced into 
SATANS-II by using the imperfection deviations along discrete 
circumferences to determine the associated meridional and 
circumferential imperfection rotations at the mesh points, 
performing a Fourier expansion on these rotations and then 
use of the Fourier coefficients. In cases where the impor- 
tant shell response harmonics and the harmonics required for 
imperfection representation exceed the limits of SATANS-II, 


a selection must be made. 
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Instructions for the use of both SATANS-I and SATANS-II 
are contained in Appendix C and complete listings of both 
programs are given in other appendices. SATANS-I requires 
210 k-bytes (53k words) of core on an IBM-360/67 computer 
using ‘single precision,' while SATANS-II requires 426 
k-bytes (107k words) of core. Imperfect shell analyses 
performed by SATANS-II in this report used a 35 station 
meridian and 25 circumferential modes. Execution times 
were an average of 60 minutes, although execution times 
of between 15 and 30 minutes can be expected for analyses 


using 15 or less modes. 


eG. IMPERFECT SHELLS 

SATANS-II has demonstrated an ability to perform a 
numerical structural analysis on the arbitrarily imperfect 
shell of revolution. Imperfections of over three shell wall 
thicknesses in depth have been handled. Increasing solution 
difficulty relates directly to increasing imperfection 
complexity and increasing shell instability. In those cases 
where extremely thin-walled shells result in perfect (imper- 
fection-free) buckling loads from SATANS that are higher 
than analytically determined buckling loads, a procedure has 
been devised where the percent reduction due to the intro- 
duction of full imperfections can be applied to the analyti- 
cally determined buckling load, and an imperfect shell 
buckling load predicted. These predictions have been shown 
to be very close to the actual imperfect shell buckling 


loads for the two specific axially loaded cylinders considered. 
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ADDENDUM 
AeoOULGITAL COMPUTER STUDY OF THE EFFECTS OF VARIOUS ARTIFICIAL 
BOUNDARY CONDITIONS AT THE POLE OF A SHALLOW SPHERICAL SHELL 


I. INTRODUCTION 


The design of many shell of revolution structures re- 
quires continuity of the structure at one or both poles. 
At such a pole, the radius becomes zero and special 
conditions must be imposed on the deflections and the stress 
resultants to assure finite strains and equilibrium. Asa 
consequence of the difficulty of incorporating these 
conditions into a structural analysis computer program, 
Many researchers model such a structure as having a small 
hole or a small rigid plug at the pole. This addendum 
investigates the applicability of the various artificial 
pole conditions by examining their effect on the stability 
of a clamped spherical cap subjected to axisymmetric and 
asymmetric pressure loads. 

_Various polar artificialities are discussed and the 
boundary equations associated with their use are given, as 
are those for the continuous pole. A determination of the 
applicability of the various poles is made by comparisons 
of their effect on the shell buckling loads and by comparing 
the errors in the normal displacement and meridional mem- 


brane force distributions incurred by their use. 
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The computer program SATANS, developed by Ball Rete 3S/ 5 
was used to conduct this study and selected results of the 
SATANS analyses are checked for validity with the computer 
programs BOSOR /Ref. 15/ and STAGS /Ref. 15/7. SATANS has a 
special routine that incorporates the proper conditions to 
be enforced at a pole. BOSOR also has this special ability 
but STAGS is not yet so Seesaeeae 

SATANS, a one-dimensional finite difference equilibrium 
formulation based on Sanders' equations, was run on an 
IBM-360/67 computer under OS/MVT. BOSOR, a one-dimensional 
finite difference energy formulation based on the Donnell 
equations, was run on a UNIVAC-1108 computer under EXEC-VIII, 
while STAGS, a two-dimensional finite difference energy 


formulation was run on a CDC-6600 computer under COMPASS-12. 
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II. STATEMENT OF THE PROBLEM 


A. GEOMETRY AND LOADING 
The structure investigated 1s the clamped shallow 
spherical cap whose geometry can be specified by the single 


nondimensional parameter i, where 


d = 2(3(1-v2))*(H/n)? (AD-1) 


and where H is the rise of the shell and h is its thickness. 
Bor this study, Poisson's ratio was taken as 0.3, Young's 
maemlus as 30 x 10° psi and h as 1.0 inches. Shells of 
X = 6 and A = 18 were investigated however only the results 
of the } = 6 shell are given in detail. 

The classical buckling pressure of a complete sphere is 


given as Gor where 
bg 2 2,73 
Geo 2B / RL (3l=V") } (AD-2) 


and where Rg rs the wadiusroft curvature... For this study, 
two loading conditions were used. The cap was subjected to 
a uniform normal pressure (symmetrical load) and a uniform 
normal pressure over half the cap (asymmetrical load). 
These loadings are shown in Figure 20. In each case, the 
magnitude of the applied pressure is q and the buckling 


results are presented as the ratio q/do- 
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A=6 SHALLOW SPHERICAL CAP 
GEOMETRY AND LOADING 


FIGURE 20 
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B. POLE CONDITIONS 

A shell is considered to have a pole if the structure is 
continuous through its axis of revolution. At a pole, the 
radius is zero and special conditions must be imposed on 
deflections u, v and w, and on the meridional moment Meee 
assure finite strains and equilibrium. The computer program 
SATANS is equipped with the capability to impose these 
special boundary conditions. A full derivation of the 
required conditions is given in Appendix D of Ref. 3. The 
results of that derivation are presented below. 


For the zeroeth Fourier harmonic (n=0): 


u=ve= Mog = t 6 = 0 
(AD-3) 
m = mM 
s6 S 
t= ta b(1-v) (u*+wowt8 )-t, 
For the first Fourier harmonic (n=l): 
(utv) = w= u* = 0 
Soe tsa. oe (AD-4) 
t 7 t6 = 56 ae 


When a double sign (+) appears, the upper sign applies for 


@n initial pole and the lower sign applies for a final pole. 


94 


















A he 

amo tT TeMoS a 
“22 a =f Sidy = 97Sn es he abbtene> et If 

A soliviovs? 29 Sine sa. Cree cee 
if seum goals phos Cetenge (tease 
noLlbtion sh? qt bes 7s aie oe «ol ee 
nanee ed? .edeteiives Bas om: >a abt 
Se > ongeaeamel, aot. 3d meee. eae ate heqgiupe " 
) o AO) Jy 2 sae, £1 me on SRE SOROD ‘iehaueds 
> | ulimegua: th Gevip "2 nnde3 £bm0% o 
Pa. od“ besas samy 296 aaltertrahs 4ads ais a 
;10°n) sia sof tcat Agate * o 


| 2 
| 
| 


ct BE enact = eg oh - 
. eS ay" \ : 
wd eT he o) (eden © - * 


« 
it=) Sic welt jestiy 


oe =o se. €.4ere; 
o = gX © gt gt 


7 Ss oe* Fae. . _ _ 
; , ~ 
a 


ign) mele Bs 
(2 # 30° Gotta te abel 


<n de ; 
7. 





For the second harmonic (n=2): 


eo A are a So 0 
pier: tate 
8 8 
i : (AD-5) 
My = rM, = +Mg, 
cae Dil Sv nia crs Bs) 
For all higher harmonics (n>2) 
u=v=we= 0 
A M = Meo = 0 (AD-6) 
Cae Cae Scion 7° 


Bore LOtations at the pole, for: the first harmonic (n=l): 


aa -w°+wu 
b9 = Fos (AD-7) 
oy Sey eu] 
and for all other harmonics, 
Cee = Gin) 10 
Sas (AD-8) 
oO) = o4(2v- = nu) 


Numerous artificial boundary conditions can be substi- 


Euted for the correct conditions at the pole. For example, 
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a small hole can be placed at the pole and can exist with 
the boundary completely free, or have conditions imposed 
such that the meridional rotation b, and deflection u are 
restrained. A meridional line load N, of magnitude suffi- 
cient to replace that 'lost" due to the hole can be applied. 
Another artificial pole condition is that best described as 
a small rigid plug. Here, both the meridional as well as 
the circumferential deflection is restrained along with the 
meridional rotation. A normal line load -P whose total 
equals that loading 'lost' due to the hole or plug can be 
applied at the boundary of the hole or plug. 

Preliminary numerical studies showed that the line load- 
ings effected a noticeable change in the solution only in 
the immediate area of the boundary, hence, the line loads 
were not considered further. Preliminary studies also 
showed that the cap behavior exhibited only minor changes 
when the small hole with a partially restrained edge was 
Substituted for that with a totally free edge. This 
investigation, therefore, concentrated on two artificial 
boundary conditions at a pole and on the continuous pole. 
The artificial pole conditions used and their associated 
mathematical expressions are: 


1. The totally free hole: 


eae tee none) 00, MA f= 0 (AD-9) 


Une) Oe Vier="" 0) <> OR 01> %. = 0 (AD-10) 
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©. TYPES OF ANALYSIS 

Three types of buckling analyses were conducted on the 
cap under each of the three pole conditions and these were: 

1. The nonlinear axisymmetric analysis of the symmet- 
rically loaded clamped spherical cap. 

2. The bifurcation buckling analysis of the symmetri- 
cally loaded clamped spherical cap with nonlinear prebuckling 
axisymmetric displacements. 

3. The nonlinear collapse analysis of the asymmetrically 
loaded clamped spherical cap. 

In the artificial pole cases, hole and plug radii were 
eaken as 0.17%, 1.70% and 3.30% of the shell meridian. The 
SATANS analyses used a 100 station meridian for the axisym- 
metric cases, a 66 station meridian with three harmonics for 
the bifurcation buckling cases and a 40 station meridian 


with five harmonics for the nonsymmetric collapse cases. 
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III. DISCUSSION 


The results of the SATANS analyses are presented in 
Table XI for all cases. Also presented are the normal 
deflections of a point located at 47.5% of the meridian 


which is near the maximum displacement for each case. 


A. THE A = 6 SHALLOW CAP 
1. The Axisymmetric Analysis 

In Table XI, the buckling load is seen as be basi- 
cally unaffected by the altered pole conditions. The normal 
displacement, however, is seen to show a marked sensitivity 
Eemene pole conditions. Figure 21 shows this sensitivity. 
In the figure, the upper curve is the deflection shape of 
Enerpcap at incipient buckling with a Cone inuous pole, and 
the lower curves represent the percent error induced in 
the normal displacement at incipient buckling by the intro- 
auetvon Of artificial poles. The 0.17% plug and the 0.173 
hole behaved identically. Holes, in general, produced an 
Beer ace 20% error in the normal displacement and this error 
Semisted Over the entire shell. The plugs, preserving the 
continuity of the structure, produced an average 10% error, 
but again, this error extended over the entire shell. 
Figure 22 shows the behavior of the meridional force Ng for 
Saeecontinuous pole case, and the N, error due to the 
eerie rcial pole conditions, also at incipient buckling. 


Large errors existed in all cases over the first third of 
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TABLE XI 


SHALLOW CAP (A=6) BUCKLING LOADS AND DISPLACEMENTS AT A POINT 





ANALYSIS POLE MERIDIONAL NORMAL 
TYPE CONDITION | STATIONS / DEFLECTION 
HARMONICS AT (47.53 
MERIDIAN 
0.9950* - 
Pole LL. 06150 1237 210 
1.0200 123200 
E52 TSUN SS LC) a a a a all a 
Buckling -0.173% 0.9720 1.1600 
- Hole-1.703% 0.9690 1.0560 
Symmetric -3.30% 0.9620 1.0290 
Load | wr rrp 
-0.173 0.9720 2570 
Plug-1.70% 0.9770 1. 1630 
-3.303% 0.9900 1.22910 
O27750* - 
Pole O273 5) 0.6386 
0.37752 0.6188 
Bifurcation §$| -o- rrr re pr rr rrr rer rrr rete 
Buckling -0.173% O37 7139 0.6356 
- Hole-1.703% 0.27635 0.6562 
Symmetric -3.30% 0.27392 0.6264 
Load | rrr re pr rrr pr rr rrr per 
-0.17% On 7 TS 05-6356 
Plug-1.70% 0.7726 0.6438 
-3.303 On7 713 0.6513 
OS 7 OOF - 
Pole 0.5391 0.9830 
BaUtapSC 9 em a a a a a a = === 
Analysis -0.173% 0353/63 0730 
- Hole-1.703 0.5389 11,0060 
Asymmetric -3.30% 0.5962 1 AAO 
Load || rr rr re rrr te pe 
-0.173 0335973 0.9478 
Plug-1.70% 0.6290 0.8100 
-3.303 027576 0.7248 


seeHuang /Ret. 16/ 


* Famili and Archer /Ref. 177 
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the meridian and were an average of 10% or less for the 
remainder of the meridian. 
2.° Biturcation Buckling Analysis 

In this analysis, the shell buckling load appears to 
beminsensitive to artificial pole conditions for the case of 
the plugs, and shows some sensitivity in the case of the 
hetes. This is surprising since bifurcation buckling of the 
Xk = 6 cap occurs in the Cos(26) mode. While equations (AD-8) 
show that ney FOtakLons “exist at whe pole) inethis: mode, rota-— 
erons certainly do exist at 1.7% and 3.3% meridian due to 
the buckling mode. The equations for the rigid plug's 
boundary conditions (equations (AD-10) ) specifically pre- 
vent this rotation, while those for the hole (equations 
(AD-9) ) do not. 

Examination of Figures 23 and 24 show that the dis- 
placement and force distributions of the plugs at incipient 
buckling to be less erroneous, again, than the poles. In 


‘ Figure 23, an average 10% displacement error occurs over the 


oo 


entrce Shell for the plugs, while holes of 1.7% and 3.3 
meridian cause large errors over much of the shell. This 
trend is repeated in Figure 24, although the errors are not 
as large. From this analysis, the indications are that the 
Metsenole is a poor artificial pole and a 3.3% hole is 
Mnacceptable as it causes a 4.6% error in buckling load and 
Over 40% error in some regions in both normal displacement 
and meridional force distributions. Larger hole sizes would 


apparently increase these errors. 
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3. Nonsymmetric Collapse Analysis 


The results for the collapse analysis are also pre- 
sented in Table XI, and in Figures 25 and 26 for the 6 = 0° 
meridian, which is symmetrically located under the loaded 
Begion. In Table XI, it is shown that the buckling loads of 
the cap with rigid plugs are quite erroneous (10.8%, 16.6% 
and 41.4% for the 0.17%, 1.7% and 3.3% plugs, respectively) 
while those for the small holes are quite reasonable. [In 
Figures 25 and 26, it is seen that the small holes cause an 
average error of 10% in both the normal displacement and 
meridional force distributions at incipient buckling over 
the entire shell, while the errors caused by the rigid 
plugs is large. 

In an asymmetrically loaded shell, large rotations 
will exist at the pole and the restrictions imposed by the 
rigid plug appear to make it an unacceptable artificial pole 
for use in this case. 

Ragures 25 and 26 also indicate that a 1.7% hole sis 
the best artificial pole for this case. Assuming such a 
"best' hole size exists for another given shell or loading 
condition, an investigation similar to that presented here 
would have to be performed to find it. This is an inherent- 
ly expensive and time consuming process and there is no 
guarantee that there is a 'best' hole size without its 
radius becoming so small that mathematical difficulties 


arise during the solution process. Consequently, the use of 
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artificial pole conditions in nonsymmetric problems must be 


approached with caution. 


B. THE A = 18 SHALLOW CAP 

A deeper shell with A = 18 was investigated in a similar 
fashion and produced similar results. The shell buckling 
load was negligibly affected (except in those cases already 
mentioned for the X = 6 shell) but displacement errors of up 
to 30% and meridional membrane force errors of an average 
10-15% were noted. In general, the X4 = 18 shell appeared 
Zee Sensitive to the type of artificial pole (except in the 
Beymmetric case), but as sensitive to any artificial pole as 


the shallower A = 6 shell. 


C. COMPARISONS WITH OTHER RESEARCHERS 

As indicated in the introduction, various SATANS results 
were checked for validity with BOSOR and STAGS. These com- 
parisons are presented with others in Table XII and very 


favorable agreement is seen to exist. 
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EV CONCLUSIONS 


tion buckling cases Livthe: radiinc Of the artificiality as 
kept at or below 3% of the meridian, In the case Of the 


asymmetrically loaded Shell, this insensitivity Only applies 


2. The displacement and force error distributions are 
Sensitive to the size and type of artificial Pole and are 
usually greater than 10% with =Omer Cases Of over 30%, and 
“they exist Over the entire Shell We Aaa. 

ae) Small UG UG Plugs appear to be the best polar models 
fOr symmetrically loaded Shells, Apparently their 
continuity-preserying characteristics greatly outweigh their 
Bo 2-tenal-restricting Characteristics, 

ae) Small holes appear to be the best polar models for 
asymmetrically loaded Shells, Here, the rotational restric— 
tions of the rigid Plug make its use unacceptable, Hor “this 
Be cha loading, an investigation will be required to deter- 


Mine the Optimum hole Size, 
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5. In general, the greater the displacement, or the 
higher the membrane force that would exist at a pole, the 
greater the error introduced by the use of an artificial 


pole. 
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APPENDIX A 


NONLINEAR TERMS 


The nonlinear terms of equations (15) are given below. 
Terms involving the initial imperfection rotation terms 


along have been omitted as they cause no straining. 
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APPENDIX B 
EIGENVALUE PROBLEM MATRIX ELEMENT DEFINITIONS 


The elements of the /A / and /B / matrices used in the 


eigenvalue solution are defined below. 


Bey. Ga y= E97 
A> Fi 22° 


Ee = FO 


Boy = Fy, 0 

Ay. = Ses Oe +1GS 5, = (G5, Gi 5/6, 4 (Be1) 

i (Gg BO" = Ges, - BO 1 /G. 

uaa he 

Boo = Gao 7 E30 ra Gin Gey 7 BawO 1/7 Gia 

4 Cay 7 Hg Os om (G5, i Be (OG) (G5 a B07 Ga, 
where 

Q = mt/L 


and where 


By, = n7/4 

B,, = no/4 

By, = 0 

Be, = n/4 

Ba = O° /4 (B-2) 
Be = 0 

Iyer ao 

B,, = 0 

Ba. = 0° 


LES 





and where all elements of the 4 x 4 matrices / E /, / F / 


and Ge. are zero except for the following, defined as 


Ei 
E22 


E23 


and where 


b 


k(1l-v)b + 9(1-v)d/(8R;) 


3 (1-v) dn/(2R¢) 
E23 
2(1-v)dn*/Ré 
1 


(5! 


4s(1-v) bn/Ry - 3 (1-v) dn/(8R6) 


bv/Rg - e(i=vy nn" 7/2, 


=e se 


= Bayes 


Sat) en / B. - (1-v) dn*/(8RQ) 
son / Ra - (1-v) (1+v)dn?/R, 


-bn/R5 - (1-v) (1+v)dn?/R, 


-vn/R, 

G23 

seer ya 
vn /Ro 
vdn/Ré 
vdn*/Ré 
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APPENDIX C 


INPUT DATA GUIDE FOR SATANS-I AND SATANS-II 


Card Columns Format Item Example Meaning 


uu i—72 
2 as) 
2 6-10 
2 Ase) 
2 6=20 
2 DZ, 5 
2 26-310 
2 gi=35 
2 36-40 
2 41-45 


18A4 


x5 


£5 


15 


1S) 


TS) 


LS 


1&5) 


TITLE - Enter any 72 characters 
NO ib The problem number, 
0<NO<10000. 
$DYNAMC F Forsa Static analysis, 
set SDYNAMC = F. 
iT For a dynamic analysis, 
set S$DYNAMC = T. 
IMODE 0 Por no modalToutput data. 
1 For modal output data. 
NDIMEN 0 Dimensional output data. 
1 Nondimensional output. 
NTHMAX 8 Summed solution will be 
printed at NTHMAX merid- 
ians, 0<=NTHMAX<=36. 
IFREQ 3 Solution will be printed 
ate the, farst) station, 
every subsequent IFREQ 
station and the last sta- 
tion, 0<IFREQ<=KMAX. 
EPRENE 2 Every IPRINT converged 
solution will be print- 
ed. 
IBCINL =]. Lf the shell. has: a pole 
ate the: first, station. 
0 If the shell has no pole 
atothe first station. 
IBCFNL -1 his the! shell has: ia, pole 
at the» last: station. 
0 Try the shell has no: pole 


atuchenwlast statvon:. 


dA) 


RE ots Tr 


“ay 





Card Columns Format Item Example Meaning 


2 46-50 I5 KMAX 3.5 Number of meridional sta- 
tions. 
Note: KMAX<201 for SATANS 
-I without plots and 
KMAX<101 for SATANS-I 
with plots or for SATANS 
=—[T. 


2 51=55 LS MNMAX 7 Number of series coeffi- 
cients used to describe 
Ehev initial conditions, 
pressure and thermal 
loads (and initial imper- 
fections if using SATANS 
-II). 

Note: For SATANS=—1: 
0<MAXM<11 and 
KMAX*MAXM<201. 

For SATANS-II: 
0<MAXM<26 and 


KMAX*MAXM<951. 
2 62-65 15 LSMAX al For a linear analysis. 
99 Use many load steps for 
a nonlinear static analy- 
sis. 
750 For a dynamic analysis, 


LSMAX is the number of 
time increments, where 
LSMAX=T,,,,,/AT. 


2 66-70 a5 LCHMAX 3 The number of load step 
size reductions ina 
static analysis, recom- 


mended range = 2-4. 
0 For a dymanic analysis. 
2. 71-75 L5 ITRMAX ul For a linear analysis. 
30 The number of iterations 


at a load or time step. 
For a nonlinear analysis, 
suggested range = 10-30, 
up. to, 50 ‘for special 
cases. 


Z 76-80 15 LG 0 initivalicconditizons. Set 
to; 0) for vaistatic analy= 
sis, or for a’ dynamic 
analysis where the shell 
isSsat inestaat.T=0'. 

il For a dynamic analysis 
Wath inatial conditions. 


118 





Card Columns Format Item Example 
3 P=t2 “EL2.3 «NU 0.3 

3 2-24 £12.3 SIGO 1000.0 
1.0 

3 24-36 E12.3 ELAST - 3E8 
IO) 

3 37-48 E12.3 TKN ~4E-2 
1.0 

3 a29=60 Hi2.s3° CHAR 8.16 
1.0 

3 61-72 E12.3 TEEO 0.0 
~996E-5 

4 Jel2 E12.3 DELOAD 0.2 
~1823E~-6 

4 13824 E12.3 EPS 0.01 


I ee ee me ee 


Meaning 
Powsson'!S ratLoy.\V. 


Reference stress level, 
Oo. 

he the: inputadata»is 
dimensional. 


Reference modulus of 
elasticity, Eo. 

If the input data is 
dimensional. 


Reference thickness, h 
If the input data is 
dimensional. 


Characteristic shell 
dimension, a. 

LG the anputiudata’ is 
dimensional. 

If a static analysis. 


Reference, time, TQ. 


For a static analysis, 
DELOAD is the load incre- 
ment. It remains un- 
changed until the solu- 
tion fails to converge 

in ITRMAX iterations, 
when it is reduced by a 
factor of five. A maxi- 


mum of LCHMAX 
ductions will 
For a dynamic 
DELOAD is the 


such re- 
occur. 

analysis, 
nondimen- 


sional time increment. 


The convergence criterion 
recommended range of 
0.0L<EPS<0.00l1. 


Card 4A is only included for a SATANS-II analysis. 


4A i= 5 t3 JUMP 1 


2 


Lg 


For an analysis using 
Single series expansions. 
For an analysis using 
double series expansions. 
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Card Column Format Item Example Meaning 


4A 5-10 15 MPERFS 0 An analysis without im- 
perfections. 
Ab An analysis with imperfec- 
tions. 


Note: if JUMP=2, MPERFS 
may be 0 or 1. 

if JUMP=1, MPERFS 
must be 0 

if MPERFS=1, JUMP 
must be 2: 


Include as many cards 5 as necessary to specify NTHMAX 
meridians. If NTHMAX equals 0, omit card 5. 


5 L=7 2, 6E12.3 LOSO A list of circumferential 
coordinates 6, in degrees 
and tenths, where the 
solution printout is de- 
Sired. The list must have 
NTHMAX entries. 


fhe TBCINL=-1, omit cards 6 through 14. If IBCFNL=-l, omit 
eaeais 15 through 23. Cards 6 through 23 describe the 
boundary conditions at the first, and then at the last sta- 
tion. The boundary conditions exist on the total variables, 
not on the individual harmonics. Loadings applied through 
specification of boundary conditions are taken in the zeroeth 
harmonic (n=0) only, as the column matrix {%} is set to zero 
for harmonics greater than zero. The boundary conditions are 
Gemensional. The format of cards 6 through 23 is 4E16.8. 





meme, 25 Card 7,16 Card 8,17 Card 9, 18 
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Gard 24 is: 
ie ineluded for ia SATANS-L static analysis. 
Ze) ineluded but blank for a SATANS-1 dynamic: analysis. 
Ss. Omitted for a SATANS-II analysis. 


Card Column Format Item Example Meaning 
24 "2 L2 SPLOTS F Indicates plots are NOT 
desired. 
T Indicates plots are 
desired. 


24 3- 4 L2 SMODAL F Indicates plots are for 
summed solutions only 
dy Indicates plots are for 


modal solutions only. 


For the remainder of card 24 entries, 0 indicates that no 
plots are desired for the particular item, and 1 indicates 
jnat\| they are desired. All graphs, are plotted as the in- 
dicated item versus the station number. If a complete plot 
is desired, insure IFREQ=l. 


24 >=) 6 ¥2 IRADII ut Plot the radii as computed 
by subroutine GEOM. 

24 7- 8 12 TGAMMA ak Plot p°/p as computed by 
subroutine GEOM. 

24 9-10 1e2 IOMEGS 1 Plot w_ as computed by 
subrouEine GEOM. 

24 P=]. 2 Ta IOMEGT it Plot w, as computed by 
subroutine GEOM. 

24 13-14 m2 IDEOMS 1 Plot wo as computed by 
subroutine GEOM. 

24 15-16 in EBGEEE Ie Plot the stifiness das 


computed by subroutine 
BDB. 
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Card Column Format Item 


24 
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24 


24 


24 


24 


24 


24 


24 


24 


24 


24 


24 


24 


24 


IL Ms: 


92,0 


Zi 2 


23-24 


25-26 


20-216 


29-30 
S1—32 


33—34 


55=36 


37-318 


39-40 


41-42 


43-44 


45-46 


47-48 


LZ 


I2 


i2 


EZ 
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T2 
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TDS TLE 


EBBSTE 


TDDSTE 


LPR 


=PS 


EP? 


dB 
IMT 


EDT TE. 


IDMT 


INS 


INTH 


INSTH 


IQS 


IMS 


IMTH 


Example 


i 


Vz 


Meaning 


Plot the stiffness d as 
computed by the subroutine 
BDB. 


Plot the stiffness db/dZ 
as computed by subroutine 
BDB. 


Plot the stiffness dd/dzt 
as computed by subroutine 
BDB. 


Plot the normal component 
of the pressure load. 


Plot the meridional com- 
ponent of the pressure 
load. 


Plot the circumferential 
component of the pressure 
Load: 


Plot the thermal load. 
Plot the thermal moment. 


Plot d/dt of the thermal 
load. 


Plot-d/de-of the thermal 
moment. 


Plot the meridional mem- 
brane force distribution: 


Plot the circumferential 
membrane force distribu- 
Eron 


Plot the meridio- 
circumferential membrane 
force \diSEri butELon. 


Plot the transverse force 
GiStri bution. 


Plot the meridional moment 
avstrrbueion. 


Pitot the carcumferential 
moment distribution. 
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Card Column Format 


24 
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24 


24 


49-50 


p52 


53—94 


55=516 


Sens: 


59-60 


61-62 


I2 


1874 
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12 


T2 


TZ 


E2 


IU 


IV 


IW 


EPHTS 


TPHET 


EPHE 


Meaning 

Plot the meridio- 
circumferential moment 
dustrribut ion. 


Plot the meridional dis- 
placement distribution. 


Plot the circumferential 
displacement distribution. 


Plot the normal displace- 
ment distribution. 


Plot the meridional rota- 
Elon distributzon. 


Plot the circumferential 


YTOtation distrpbution. 


Plot the meridio- 
circumferential rotation 
aLSeributron. 


Insert imperfection data here for a SATANS-II analysis with 


imperfections. 


Insure format of the imperfection data is 


compatible with that specified in the user-written sub- 
routine IMPERF. 


ee re ee me me ee ee ee me es ee ee ee me ee ee ee ee ee we ew we 


ZS 


Indicates this is the 
onlyenun. 

Indicates another run is 
to be made. Add another 
complete set of data 
Cards carter his’ card. is 
IRNAGN=1. 
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LISTING OF SATANS-II MODIFIED SUBROUTINES AND COMMON STATEMENTS 
UBROUTINES WHERE THEY OCCUR. 
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